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ABSTRACT 

The Vela pulsar is the brightest persistent source in the GeV sky and thus is the traditional first target for 
new 7-ray observatories. We report here on initial Fermi Large Area Telescope observations during verification 
phase pointed exposure and early sky survey scanning. We have used the Vela signal to verify Fermi timing 
and angular resolution. The high quality pulse profile, with some 32,400 pulsed photons si E > 0.03 GeV, 
shows new features, including pulse structure as fine as 0.3 ms and a distinct third peak, which shifts in phase 
with energy. We examine the high energy behavior of the pulsed emission; initial spectra suggest a phase- 
averaged power law index of F = 1.5lj;QQ4 with an exponential cut-off at Ec = 2.9 ± 0.1 GeV. Spectral 

fits with generalized cut-offs of the form er^^/^"^'' require 6 < 1, which is inconsistent with magnetic pair 
attenuation, and thus favor outer magnetosphere emission models. Finally, we report on upper limits to any 
unpulsed component, as might be associated with a surrounding synchrotron wind nebula (PWN). 
Subject headings: gammaiays: general; pulsars: individual: PSR B0833— 45 
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1. INTRODUCTION 

Radio pulsations at P=89ms from PSR B0833-45 (=PSR 
J0835— 4510) in the Vela supernova remnant were discovered 
by Large, Vaughan & Mills (1968). This pulsar is bright 
(51.4 GHz ~ 1.5 Jy), young (characteristic age Tc = P/2P = 
llkyr), and energetic (Esd = 6.9 x 10^^/45 ergs/s, for a 
neutron star moment of inertia 10^^/45 gcm^). It is embed- 
ded in a flat spectrum radio synchrotron nebula (see Frail et 
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al. 1997) and is surrounded by a bright X-ray wind nebula dis- 
playing remarkable toroidal symmetry (Helfand et al. 2001; 
Pavlov et al. 2003). A VLBI parallax measurement provides 
a well-determined distance of £> = 287^ J^pc (Dodson et al. 
2003), improving on earlier optical HST measurements (Car- 
aveo et al. 2001). This proximity means that the spindown 
energy flux at Earth, Esd/^t^<P, is second only to that of the 
Crab pulsar and ensures that Vela is among the most intensely 
studied neutron stars, particularly at high energies. 

Pulsed 7-ray emission from Vela was first detected during 
the SAS-2 mission (Thompson et al. 1975); it is, in fact, the 
brightest persistent source of celestial 7-rays and has been the 
proving ground of GeV observatories ever since. The ba- 
sic source properties discovered by SAS-2 were elaborated 
using observations by COS-B (Kanbach et al. 1980), and 
CG/?0/EGRET (Kanbach et al. 1994): the source is ~100% 
pulsed with no convincing evidence for year-scale variations, 
the two 7-ray peaks are separated by 0.42 in pulsar phase and 
the first peak lags the radio peak by 0. 12 in phase. The source 
spectrum is hard with an average photon index F = 1.7 with 
evidence of a cut-off above 2-A GeV. The most detailed study, 
to date, of the 7-ray pulsations is that of Fierro et al. (1998), 
which produced a high signal-to-noise pulse profile and evi- 
dence for phase-resolved spectral variations. As the present 
work was being prepared for submission, the AGILE team 
reported their initial results on 7-ray pulsars, including Vela 
(Pelli/oni et al. 2008); the results, given the more liinited 
count rates and energy resolution, are broadly consistent with 
our conclusions described below. 

With the successful launch of the Fermi Gamma-ray Space 
Telescope, formerly GLAST, observatory on June 11, 2008, 
we have a new opportunity to examine the high energy be- 
havior of the Vela pulsar and to study this archetype of the 
young, energetic pulsars in detail. During Launch and Early 
Operations phase (L&EO) the Fermi mission targeted the Vela 
pulsar for a number of pointed observations, in addition to 
coverage during initial tests of the sky survey mode. In the 
latter mode the instrument axis is offset North and South of 
the zenith during alternate orbits to provide near-uniform sky 
coverage every three hours. One main purpose of these early 
observations was to tune the Large Area Telescope (LAT) per- 
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formance on celestial 7-ray sources. However, the initial re- 
sults on Vela itself are of interest, including new high energy 
features in the pulse profile, an improved measurement of the 
high energy cut-off in the pulsar spectrum, and a search for 
associated pulsar wind nebula (PWN) emission at GeV ener- 
gies. 

2. RADIO TIMING OBSERVATIONS 

The Vela pulsar is young and exhibits substantial timing ir- 
regularities. This means that the optimal use of the Vela 7-ray 
photons, which arrive at an orbit-averaged rate of one every 
4 minutes during LAT sky survey observations, requires a si- 
multaneous radio ephemeris. The radio ephemeris is obtained 
using observations made with the 64-m Parkes radio telescope 
as part of the overall program for pulsar monitoring in support 
of the Fermi mission (Smith et al. 2008). A total of 27 times- 
of-arrival (TOAs) were obtained at a frequency of 1.4 GHz 
over the period of the LAT observations. The median error in 
each TOA is 1 .7 /xs. Fits to the TOAs were made with the pul- 
sar's rotational frequency and frequency derivative as free pa- 
rameters and the ephemeris created from these data maintains 
phase with an rms residual of 90 /is, or 10^'^ of pulsar phase, 
throughout the LAT observations. No large timing irregular- 
ities or glitches were detected. An accurate determination of 
the pulsar's dispersion measure of 67.95 ± 0.03 cm~^pc was 
obtained by measuring the delay in the TOAs of the pulse 
across the 256 MHz bandpass of the 1.4 GHz receiver. This 
allows extrapolation of the radio ephemeris to infinite fre- 
quency with an error of ~60 /zs. Photon arrival times were 
referred to the Solar-System barycenter and pulse phases were 
assigned using the standard pulsar timing software TEMP02 
(Hobbs et al. 2006). 

We thus can assign a pulsar phase to the 7-ray photons ref- 
erenced to the radio with high confidence. Gamma-ray events 
recorded with the LAT have timestamps that derive from a 
GPS clock on the Fermi satellite. Pre-launch ground timing 
measurements using cosmic muons seen coincidentally in the 
LAT and in an independent counting system demonstrated 
that the LAT measures event times relative to the spacecraft 
clock with a resolution of better than 300 ns. On orbit, satel- 
lite telemetry indicates a comparable resolution. The contri- 
bution to the barycenter times from uncertainty in the LAT's 
position is negligible. It is also useful to confirm the absolute 
timing using celestial signals; as of this writing the six clas- 
sical pulsars detected by EGRET have provided good quality 
LAT pulse profiles. The phasing of the radio and 7-ray peaks 
is consistent with the best EGRET results, ensuring absolute 
timing relative to the radio time system better than ~ 1ms. 

3. GAMMA-RAY OBSERVATIONS 

The LAT instrument on Fermi is described by Atwood et 
al. (2008) and references therein. The LAT is an electron- 
positron pair production telescope featuring solid state silicon 
trackers sensitive to photons from < 30 MeV to > 300 GeV. 
It has a large ~2.4sr field of view, and compared to earlier 
7-ray missions, has a large effective area (~ 8000 cm^ on 
axis), improved resolution (~ 0.5°, 68% containment radius 
at 1 GeV for events collected in the 'front' section with thin 
radiator foils) and small dead time (^ 25 /is/event). We report 
here on the LAT's initial observations of the Vela pulsar, using 
data collected during 35 days of on-orbit verification tests and 
the initial ^ 40 days of the on-going first year sky survey. 
These data already suffice to illustrate the power of the LAT 
for astronomical observations and, indeed, show several new 



features in the radiation of this well-known 7-ray pulsar. 

Table 1 contains a journal of the Vela coverage, along with 
numbers of photons with measured energy above 0.03 GeV 
recorded within a radius of 5° of the pulsar, including back- 
ground photons. For both pointed and survey observations 
we exclude events within 8° of the Earth's limb to minimize 
contamination by 'Earth albedo' photons. 

During the L&EO period, the instrument configuration was 
being tuned for optimum performance. Accordingly, the 
knowledge of the energy scale and effective area are more 
limited than for routine operations. We discuss how these data 
verify the LAT photon selection, effective area, timing, pho- 
ton energy measurement and the variation of the point spread 
function (PSF) with energy. These results may be of use to re- 
searchers seeking to predict Fermi capabilities during longer 
exposures of fainter sources, in both sky survey and pointed 
mode. 

4. RESULTS 

4.1. Pulse Profiles 

After calibration of the LAT pointing axis on a large set of 
identified high latitude 7-ray sources, we find that the best fit 
position of the Vela point source is within 0.5' of the radio 
pulsar position with a statistical error of 0.4'(95%); the LAT 
has unprecedented accuracy for localization of bright hard 7- 
ray sources, although some systematic uncertainties remain. 

The Vela pulsar is embedded in the bright 7-ray emission 
of the Galactic plane, which is particularly strong at low en- 
ergies. Further, the LAT, like all pair production telescopes, 
has an angular resolution dominated by scattering at low pho- 
ton energies ^gs ~ ^■^°^Gev^ (^£6 below). Thus selection 
of photons from an energy-dependent region of interest (ROl) 
around Vela is important; the best selection depends on the 
desired product. Here, we seek pulse profiles with good sig- 
nal to noise over a broad energy range, so we select photons 
within an angle 9 < Max[1.6 — 3Logio(ii^Gcv), 1-3] deg of 
the pulsar position. This includes a larger fraction of the PSF 
at high energies, where the background is relatively faint. We 
use 'Diffuse' class events, those reconstructed events having 
the highest probability of being photons. 

In this energy-dependent aperture we have 32,400±242 
pulsed photons and 2780 ± 53 background photons with mea- 
sured energy > 0.03GeV. Figure 1 shows the 0.1-lOGeV 
pulse profile from this ROI with the peak of the radio pul- 
sar signal at phase 0. To highlight the fine structure we plot in 
Figure 1 the pulse profile using variable-width bins, each con- 
taining 200 counts. These counts, divided by the bin width, 
give the photon flux in each phase interval; these phase bin 
fluxes thus have a la Poisson statistical error of 7%. The 
pulse profile is normalized to 100 at the pulse peak. Three in- 
sets show structure near the first peak (PI), second peak (P2) 
and 'off-pulse' window at a finer scale. These peaks are at 
pulsar phases (j) = 0.130 ± 0.001 and = 0.562 ± 0.002, 
respectively. 

The peaks are asymmetric. In particular P2 has a slow rise 
and a fast fall. PI also has a steep outer edge with the fall 
somewhat slower. We fit each peak in Figure 1 with two 
half-Lorentzian functions (with different widths for the lead- 
ing and traihng sides) over phase intervals which avoided the 
complex bridge flux (0.11 < < 0.16 for PI, 0.53 < ^ < 
0.59 for P2). The outer edges of the two peaks had consistent 
Lorentzian half- widths of (p — 0.012±0.001. Extra structure 
at the pulse peaks made the half maximum widths somewhat 
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TABLE 1 

Early LAT observations of PSR B0833-45 



Date 


MJD 


Primary mode 


Ar(E>0.03GeV)" 


Notes 


2008 Jun30-Jul4 


54647.40-54651.36 


Sky survey 


1859 


LAT First light. 


2008 Jul 4-Jul 15 


54651.38-54662.10 


Survey+Pointings 


5172 


LAT Calibrations. 


2008 Jul 15-Jul 19 


54662.12-54666.08 


Pointings+Survey+limb FoUowing 


1170 


Pointed Obs. Tuning 


2008 Jul 19-Jul 22*'' 


54666.10-54669.14 


Pointed Vela+ 2"'' Target 


7751 


Pointed Obs. Tuning 


2008 Jul 22-Jul 24 


54669.16-54671.41 


Pointed Vela+ 2"<* Target 


3212 


Pointed Obs. Tuning 


2008 Jul 24-Jul 30* 


54671.36-54677.45 


Pointed Vela+ 2"'' Target 


7607 


Pointed Obs. Tuning 


2008 Jul 30-Aug 3 


54677.45-54681.66 


Polar study + Sky Survey Tuning 


886 


Nominal Ops. 


2008 Aug 3-Sep 14* 


54687.68-54723.91 


Sky Survey 


26680 


Nominal Ops.'^ 



"AH photons (pulsed and background) within 5° of PSR B0833-45, zenith angle z <105°, event class Pass6-Dijfuse. 
''Observations marked by asterisk are those used for spectral studies. 
'^L&EO ended on Aug. 1 1 . 




Phase 

Fig. 1 . — Vela broad band (B = 0. 1 — 10 GeV) pulse profile for all photons from an energy dependent ROI. Two pulse periods are shown. The peak of the 
radio pulse is at phase = 0. The coimt rate is shown in variable width phase bins with a constant number of counts. The dashed line shows the backgroimd 
level, as estimated from a surroimding annulus during the off-pulse phase. Insets show the pulse shape near the peaks and in the off-pulse region. 
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smaller ((f> = 0.009 ± 0.002; 0.8 ms). The falling edge of 
peak 1 has a Lorentzian half width of (p = 0.017 ± 0.0015, 
while the rising edge of P2 has a width = 0.027 ± 0.005. 
Both peaks show apparently significant structure on scales as 
small as S4> « 0.003 (~ 300 /us), but additional counts are 
required to fully probe the pulse profile at this scale. We also 
show an inset of the 'off pulse region' ((f) = 0.65 — 1.05), 
with 50 counts/bin to show fine structure. Here and in the 
main panel, the dashed line shows the estimated residual un- 
pulsed background counts in our energy-dependent ROI, as 
measured from the pulse minimum. This is in good agree- 
ment with a background level estimated from 3 — 5° around 
the pulsar. We see that, as for other non-radio Vela pulse pro- 
files, there is a fiiint tail of pulsed emission in this region, 
reaching non-detectability only near ^ = 0.8 — 1.0. We es- 
timate that this 'off pulse' window contains 235±15 pulsed 
photons (7.3 X 10"'' of the pulsed flux). 

There is dramatic energy evolution in the 7-ray pulse pro- 
file. Figure 2 shows pulse profiles from six energy bins, cov- 
ering three decades of the LAT data, drawn from our energy- 
dependent aperture. Note that the relatively large aperture 
at low energies leads to larger off-pulse flux from the sur- 
rounding background. We should also caution that, with 
the present limited accuracy in energy measurement below 
0. 1 GeV, some higher energy photons may leak into the lowest 
energy bin. The most prominent pulse feature is a decrease of 
PI relative to P2 with increasing energy; PI is not detectable 
above ~10GeV (Figure 2). This confirms a trend seen in the 
EGRET data for the Crab, Vela and Geminga pulsars - the 
second pulse dominates at the highest energies (Thompson 
2001). Interestingly, at the lower energies, below ~ 120MeV, 
the trend is reversed with PI weakening again with respect to 
P2. We do not find any statistically significant evidence for 
shifts in the phases of the narrow PI and P2 pulse components 
with energy. These structures, spanning together < 0.07 of 
the rotational phase, contain emission dominating over three 
decades of the pulsar photon output. 

The 'bridge' region between the peaks also shows appre- 
ciable structure with a trailing shoulder of PI shifting to later 
phase. At > 1 GeV it is clear that this is a distinct pulse com- 
ponent (P3) with a peak at </> ~ 0.27 in the 3-10 GeV band 
(Figure 2). Note that this peak shifts in phase by Scj) w 0.14 
between 0.2 and 15 GeV (Figure 3). 

To compare with the lower energy structure, we show at 
bottom left of Figure 2 the 8-16keV non- thermal X-ray pulse 
measured by RXTE (Harding et al. 2002). On the right we 
similarly show the 4. 1-6.5 eV NUV HST STIS/MAMA pulse 
profile of Romani et al. (2005). The 1 .4 GHz radio pulse pro- 
file, whose peak defines phase (f) = 0, runs across the bot- 
tom two panels. Both the optical/UV and the hard X-rays are 
dominated by non-thermal magnetospheric emission. In con- 
trast, the < 1 keV soft X-ray flux is dominated by thermal 
emission from the neutron star surface (e.g. Manzali, de Luca 
& Caraveo 2008), making it much more difficult to trace the 
non-thermal peaks through this intermediate energy band. 

We may compare components in the UV and hard X-ray 
pulse profiles shown with those of the non-thermal 7-ray pul- 
sations. In particular, we should note that the 7-ray PI compo- 
nent is dominant in the non-thermal X-rays. This component 
is absent in the optical/UV, but pulse profiles at these ener- 
gies have a strong peak in the bridge region at ^ ^ 0.25, well 
matched in phase to the P3 structure in the > GeV pulse pro- 
files. This possible connection between optical pulses and the 
7-ray bridge emission was already noted by Kanbach et al. 



(1980). Note that the profile for UV energies has a distinct 
sharp pulse coincident with the radio at (/> = 0. While the 
general slope of the faint 4> = 0.65 — 1.05 7-ray emission 
matches that of the UV, no sharp 7-ray pulse components are 
yet visible in this phase interval. 

4.2. The PSF and a Search for Unpulsed Emission 

Figure 1 makes it clear that the Vela extraction region is 
strongly dominated by the pulsed emission. However there 
are substantial phase windows where the pulsar is very faint. 
Figure 4 presents Fermi 7-ray images of the Vela pulsar in 
three energy bands in the LAT range separated into the on- 
{(j) = 0.05 - 0.65) and off- (0 = 0.65 - 1.05) pulse phase re- 
gions. The images (except for the low count high energy off- 
pulse frame) are smoothed by Gaussians (with FWHM = 
1.2°, 0.7° and 0.2° for the low to high energy images, respec- 
tively). The off-pulse images show the bright diagonal band 
of the Galactic diffuse emission, with appreciable structure. 
Comparison with the PSF contours (red) show that this is re- 
solved in the vicinity of Vela at the higher energies. While a 
clump of extended emission surrounds the pulsar position on 
the few degree scale, it is at present not possible to associate 
this flux with known Vela PWN structures (e.g. the Vela-X 
TeV PWN, contours in the middle energy band). 

We have subtracted these off- pulse images (scaled by 1.5 x 
to normalize exposure) from the on- pulse images and plotted 
(error bars) the profile of the resulting point source (bottom 
panels). The lines show simulations of the expected PSF from 
pre-launch calibrations ('Pass6'), computed for a source with 
a Vela-like spectrum and off-axis angles similar to the actual 
observations. The agreement is quite good. 

As a first attempt at constraining unpulsed (e.g. PWN) 
emission from Vela we have attempted to fit for a point 
source in the off-pulse phase window, fixed at the position 
of file pulsar. Using 0.1-10 GeV photons in an 8° ROI in flie 
(j) = 0.65 — 1.05 phase interval, we derive a 95% CL upper 
limit on the flux of 1.8 x 10~^7/cm^/s. After subtracting 
the estimated remnant pulsed flux in this window (0.73% of 
the phase-averaged flux) and scaling to the full pulse phase, 
this provides a limit on an unpulsed point source at the posi- 
tion of Vela of 2.8% of the E > 100 MeV pulsed emission 
count-rate. More photons, especially at higher energy, will 
be required to search for a resolved PWN correlated with the 
TeV or radio structures. 

4.3. Energy Spectrum 

To study the phase-averaged spectrum of Vela, we used the 
standard maximum-likelihood spectral estimator 'gtlike' to be 
provided with the Fermi SSC science tools. This fits a source 
model to the data, along with models for the isotropic (in- 
strumental and extragalactic) and structured Galactic back- 
grounds. The instrumental background was comparable to 
estimates from pre-launch simulation (Atwood et al. 2008). 
We used data from the observation spans indicated in Table 1, 
selecting photons with E > 0.1 GeV within 15° of Vela. Our 
basic model for the spectrum of Vela is a simple power-law 
with an exponential cut-off. 

With the large number of events collected for Vela, the sta- 
tistical errors are very small. Systematic errors for the LAT 
are still under investigation (Abdo et al. 2009). Here we adopt 
conservative estimates of the systematic uncertainty in the 
LAT effective area, derived from the on-orbit estimation of 
the photon selection efficiency as function of energy and off- 
axis angle. This estimation has been made by comparing Vela 
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Fig. 2. — The evolution of the Vela 7-ray pulse profile over three decades of energy. Each pulse profile is binned to 0.01 of pulsar phase, and dashed lines show 
the phases of the PI and P2 peaks determined from the broad band light curve. In the top right panel we label the main peaks PI, P2 and P3. In the bottom panels 
we show at left the 8-16 keV RXTE pulse profile of Harding et al. (2002) along with the radio pulse profile (in red). At lower right, the 4.1-6.5 eV HST/STIS 
NUV pulse profile of Romani et al. (2005) is shown for comparison. 

the LAT response. We conducted this exercise for a range of 
energy and off-axis angle bins, and use the precision of the 
resulting effective area estimates and the difference from the 
simulated effective area as an estimate for the uncertainty in 
the LAT effective area as a function of incident photon energy. 
This varies from < 10% near 1 GeV to as much as 20% for 
energies below 0.1 GeV and 30% for energies greater than 10 
GeV. 

The gtlike fit for 0.1 < E < 30 GeV is unbinned, and 
results in a spectrum of the form 



0.3 
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0.1 



diV 

dE 
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Fig. 3. — The evolution of the P1/P2 pulse peak ratio (open circles) and the 
P3 pulse peak phase (sohd squares) with energy. 

on-pulse event selection to event selection in the off-phase 
bins. This allows an estimate of the efficiency for a pure pho- 
ton sample, which is compared with our current simulation of 



(2.08±0.04±0.13)xl0"^.B^e"'^/^<=7/cmVs/GeV, 

(1) 

with E in GeV, T = -1.51 ± 0.01 ± 0.07 and Ea = 
2.857 ± 0.089 ± 0.17GeV. The first errors are the statisti- 
cal values for the fit parameters, while the second errors are 
our propagated systematic uncertainties. The latter are par- 
ticularly large for our spectral index parameter F because of 
the large uncertainty at present in the effective area at high 
and low energies. To explore the stability of this result, we 
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Fig. 4. — On- ((f> = 0.05 — 0.65, top row) and off- {(f> = 0.65 — 1.05, middle row) pulse phase images in three energy bands. Each panel spans 10° X 10° in 
equatorial coordinates centered on Vela (RA and Dec marked in degrees in the central panel). On-pulse data are dominated by the point source; the panels have 
a hard square root stretch to show the faint wings of the PSF. The off-pulse images (linear grey scale) are dominated by Galactic diffuse emission; some of the 
structure is due to limited count statistics. Contours of the point source at 0.5, 0.25 and 0.125 of the peak intensity are shown (red) in each frame. In the middle 
energy off-pulse image we also show the HESS Vela PWN contours (Aharonian et al. 2006; blue). Along the bottom row background-subtracted Vela PSFs (red 
eiTor bai's give the statistical eiTor in the observed counts) are compared with simulations based on pre-launch descriptions of expected LAT performance (solid 
lines) and estimates of the systematic uncertainty in the model PSF (dashed lines). 



also fit using three other spectral estimators. We used a stan- 
dard XSPEC analysis with our best model response matrices, 
a binned maximum likelihood estimator which computes the 
on-pulse photon counts in a point source weighted aperture in 
excess of off-pulse background counts ('ptlike') and a method 
which propagates the model spectrum through simulated in- 
strument response to compare with observed pulsed source 
counts. For additional tests of the stability of the spectral fits, 
the data were also analyzed separately for 'pointed' and 'sur- 
vey' observations, as well as for events in the front (thin ra- 
diator foil) and back (thick radiator foil) sections of the LAT 
detector. The fit parameters for the various data sets and anal- 
ysis techniques, do have some statistically significant varia- 
tion, but all were well within our presently-estimated system- 
atic errors, as listed in Equation (1). We believe that as our 
understanding of the LAT performance improves we should 



be able to substantially decrease these systematic uncertain- 
ties. 

Figure 5 presents the spectral power ^ along with this 
best fit model. The binned spectral points are drawn from the 
ptlike analysis and show both the statistical error flags (bars 
with caps) and the inferred systematic error flags (bars with- 
out caps); the systematic errors dominate for all energies be- 
low 7 GeV. We also plot the EGRET data points of Kanbach 
et al. (1994). These are discrepant with the LAT points, par- 
ticularly above a few GeV. Some studies have indicated that 
the EGRET response was incorrectly estimated (Baughman 
2007, Stecker et al 2008). 

We have also attempted to fit the data with a general- 
ized cut-off model of the form Exp[—{E/Ec)'^]. We find 
b = 0.88±0.04t[Jj^ so that models with a hyper-exponential 
behavior are well excluded. Taking into account the system- 
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Fig. 5. — The phase-averaged Vela spectral energy distribution (E^dNj/dE). Both statistical (capped) and systematic (uncapped) errors are shown. We 
believe that the latter are conservative; they dominate at all energies below 7 GeV. EGRET data points (diamonds, Kanbach et al. 1994) are shown for comparison. 
The curve is the best-fit power law with a simple exponential cut-off. 
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atic errors the spectrum is fully consistent with the simple ex- 
ponential 6=1 cut-off. To check the exclusion of a given 
hyper-exponential model we compute the probability of in- 
correct rejection of a given b value using the likelihood ratio 
test. For example, if only statistical errors are included 6=2 
is rejected at 16.5(t; inflating the errors to the level of our es- 
timated systematic uncertainty leads to a small, but not negli- 
gible, 0.29% chance of incorrectly rejecting the 6=2 model. 

The analyses of Kanbach et al. (1980), Grenier et al. (1988) 
and Fierro et al. (1995) already showed that there is signifi- 
cant spectral variation through the pulse, with the bridge re- 
gion harder than the peaks. Clearly, in Figure 2 we see that 
the situation is even more complicated, with three or more 
overlapping physical components contributing to the energy- 
dependent pulse profiles. Also, in Figure 5, one sees that there 
may be emission in the highest energy bins in excess of a 
simple exponentially cut-off model. This likely arises from 
the extended hard spectrum of the P2 component, and drives 
the best-fit value for 6 below 1 . Clearly the phase-averaged 
spectrum is a superposition of emission from a variety of dis- 
tinct physical regions, and we cannot expect a simple cut-off 
power-law spectral model to be an exact representation of the 
data. Reduced systematic errors and the photons of a year 
or more of LAT data will be required to explore the phase- 
resolved spectral structure. 

5. DISCUSSION 
5.1. Light Curve 

The early LAT data have already greatly clarified the 7-ray 
emission of the Vela pulsar and, even before detailed phase- 
resolved spectral measurements and model comparisons, sig- 
nificantly constrain the origin of high energy pulsar emission. 
Pulsar particle acceleration is believed to occur in the open 
zone, on the field fines above each magnetic pole which ex- 
tend through the light cyUnder. These subtend ~ 0.5 of pul- 
sar phase at high altitude, so with a total phase in the narrow 
pulses of ^ 0.06 one would expect a comparable fraction of 
the open zone, and hence fraction of the spin-down power, to 
contribute to the 7-ray emission. The narrow peaks are often 
interpreted as caustics (Morini 1983; Romani & Yadigaroglu 
1995, Dyks et al. 2004), where the combination of field re- 
tardation, aberration and time delay causes emission from the 
boundary of an acceleration zone (gap) to pile up in pulsar 
phase. A simple pulse from a 1 -D fold caustic would increase 
to a maximum as I{(f>) oc (0 — (j)o)~^^'^9[4> — 0o] as one ap- 
proaches the destruction of an image pair at phase (po or fall 
off in mirror form after image pair creation. Vela's P2 struc- 
ture with a slow rise and abrupt fall looks Hke a trailing (image 
pair destruction) caustic, albeit with a steeper increase to the 
maximum, suggesting intensity variation along the emission 
surface or a higher-order catastrophe. The shape of PI is less 
clear, but the apparent fast rise and slow fall would indicate a 
leading (image pair creation) caustic (Figure 1). 

The broad structure of the P3 interpulse region suggests 
a beam of emission from a bundle of field lines bracketed 
by the two caustic-forming zones. The similarity between 
the 7-ray P3 phase structure and that of the lower frequency 
optical/UV emission at « 0.2o suggests a common ori- 
gin - they both cover similar phase intervals and the peak 
shifts to later phase at higher energy. One appeaUng possibil- 
ity is to associate the optical/UV P3 with synchrotron emis- 
sion from relatively low altitude pair cascades produced by 
the few-GeV curvature radiation dominating the uFi, spec- 



tral energy distribution (SED). These 77 — > e^ pairs, with 
electron energy 7e ^ 2GcV/meC^ ^ 10'^-^, will produce 
synchrotron emission with characteristic peak energy e ~ 
12Si27eSin*_i keV as they radiate away their pitch angle 
~ 0.1'f_i radians in a magnetic field of 10^^ G. For Vela, 
we have B12 ~ 3(r/_R*)~^, so to produce e ~ 5 eV pho- 

— 1/3 

tons the piiirs would reside at r ^ 804'_j^ i?*, where i?* is 
the neutron star radius. If the synchrotron emission is outward 
beamed, along with the pairs, then the Compton scattered flux 
would appear at ^ 5'7pe ^ 5GeV, matching the LAT P3 
component. This suggests that the pulse of seed photons at 
IR to UV energies should have energy-dependent phase shifts 
tracking their Compton-up-scattered P3 progeny. 

5.2. Luminosity and Beaming 

While the VLBI parallax provides a well-determined dis- 
tance of _D = 287li[7pc, the pattern of the 7-ray beam and 
hence the true 7-ray luminosity are uncertain. To compute the 
true luminosity 

= 4^/nFE,obs^' (2) 

we need a correction factor /n from the observed phase- 
averaged energy flux FE,obs at the Earth line of sight (at angle 
d^E to the rotation axis) to an average over the sky. For a given 
model, the correction depends on the geometry, so we have 

fn{a,CE)= J F^{a,C,ct>)dn/{F,,p47r), (3) 

where the model for a given magnetic inclination a from the 
rotation axis gives Fj(a, the pulse profile as a func- 
tion of phase (j) seen at viewing angle C- Fexp is the ex- 
pected phase-averaged flux for the actual Earth line-of-sight 
(e- In general highly-ahgned models, such as the polar cap 
scenario, have small sky coverage and fn < 0.1, while 
outer-magnetosphere models have much larger /n (Romani 
& Yadigaroglu 1995). For Vela, the viewing angle is esti- 
mated to be C ~ 63° from the geometry of the Chandra X- 
ray (9/7.servflfory-measured X-ray torus (Ng & Romani 2008). 
The magnetic impact angle (3 = — a is constrained by the ra- 
dio polarization data. Together, these aUow us to compute the 
high energy beam shape (assuming uniform emissivity along 
gap surfaces). For the outer gap (OG) picture we compute 
/n = 1.0 for ^ = 64°. For the two-pole caustic (TPC) model 
appropriate to the slot gap picture, we find fn = 1.1 (Watters 
et al. 2008). 

The observed energy flux obtained by integrating times 
the photon number flux in Equation (1) over the range 
0.1-lOGeV is (7.87±0.33 ± 1.57)x lO'^ergs/cm^/s. For a 
distance of 0.287 kpc, the 7-ray luminosity of Vela is then 
7.8 X 10*^^/0 ergs/s. Thus the observed 7-ray flux implies an 
efficiency of r]^ = 0.011/n//45 with /n > 1 for both the OG 
and TPC models. Note that this is a substantial fraction of 
the geometrical bound on the efficiency from the peak widths 
noted above. 

5.3. Gamma-Ray Pulsar Models 

Three general classes of models have been discussed for 7- 
ray pulsars. In polar cap models (Daugherty & Harding 1994, 
1996; Stumer & Dermer 1995; Harding & Muslimov 2004), 
the particle acceleration and 7-ray production takes place in 
the open field line region within one stellar radius of the mag- 
netic pole of the neutron star. In outer gap models (Cheng, Ho 
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& Ruderman 1986; Romani 1996; Hirotani 2005), the inter- 
action region lies in the outer magnetosphere in vacuum gaps 
associated with the last open field line. Other recent mod- 
els hark back to the 'slot-gap' picture of Arons (1983), with 
the polar cap rim acceleration extending to many stellar radii 
(Harding & Muslimov 2004). 

The appreciable offset from the radio peak, coupled with 
the presence of only one radio pulse, suggests that the 7-rays 
arise at high altitude. This is also supported by the lack of 
hyper-exponential absorption in the spectrum, which is a sig- 
nature of 7-B pair attenuation of low altitude emission. The 
remaining two geometrical models can produce the general 
double peak profile of Vela (and many other 7-ray pulsars). 
In 'Outer gap' models emission starts near the 'null charge' 
Q ■ B ~ surface at radius tnc and extends toward the 
light cylinder tlc = cP/2'k. One magnetic pole dominates 
the emission in each hemisphere and the two peaks represent 
leading and trailing edges of the hollow cone of emission from 
this pole. If on the other hand emission extends well below 
the nuU charge surface toward the neutron star, both magnetic 
poles can contribute toward emission in a given hemisphere. 
This is the two-pole caustic model of Dyks & Rudak (2003) 
which might be realized in 'slot gap' acceleration models. In 
this case the leading 'image pair creation' caustic from high 
altitudes should not be visible, and the first 7-ray pulse repre- 
sents the trailing caustic from emission at r < r^c below the 
null-charge surface (the radio magnetic pole) while the second 
pulse represents a trailing caustic at higher altitudes from the 
opposite pole; both pulses should generally show a slow rise 
and steeper fall. Detailed 7-ray pulse profiles and, especially, 
phase resolved spectra can help distinguish these scenarios. 

Another geometric difference in the two scenarios is the ori- 
gin of the emission outside the main pulse. Low altitude field 
lines in the two-pole model tend to produce emission along 
all Unes of sight. In the outer gap picture high altitude field 
lines can contribute faint emission extending through the off 
pulse region. 

Further, as noted above, the energy spectrum also presents 
challenges for near-surface emission. In the polar cap mod- 
els, a sharp turnover is expected in the few to lOGeV en- 
ergy range due to attenuation of the 7-ray flux in the mag- 
netic field (Daugherty & Harding 1996). The spectral change 
at ^2.9 GeV does not appear to fit this model. We conclude 
that low altitude radiation subject to 7-B pair production can- 
not account for the bulk of the Vela 7-ray emission. Indeed 
we can use the observed cut-off to estimate a minimum emis- 
sion height as r « (f„,,„,,,.Si2/1.76GcV)2/'^p-i/7^^^ where 

€max is the unabsorbcd photon energy, P is the spin period 
and the surface field is 10^^Bi2 G (Baring 2004). Using Vela 
parameters {P = 0.089 s, B12 = 3.4) we see that our cut-off 
implies that the bulk of the emission arises from > 2.2i?*. 
Since we see pulsed photons up to emax = ITGeV, this emis- 
sion must arise at r > 3.8i?*. A similar conclusion has been 
recently made for the Crab pulsar by the MAGIC team, who 
observe pulsed E > 25 GeV emission, using the imaging air 
Cerenkov technique (Aliu et al. 2008). Finally, it should be 
noted that while typical radio pulsar emission is inferred to 
arise at relatively low altitudes, for Vela-type pulsars recent 
radio models infer emission heights of as much as IGOi?* 
(Karastergiou & Johnston 2007). All of these factors point 
to a high-altitude origin of the Vela 7-ray emission. 



6. SUMMARY 

The early LAT observations of Vela serve to show that 
the instrument is performing very well, with effective area 
comparable to expectations, good PSFs and excellent source 
localization, especially at high energies. The time-tagging 
of events is excellent and our pulse profile reconstruction is 
limited primarily by the accuracy of the radio-derived pul- 
sar ephemeris. This presently has an accuracy of ~100yus 
over our analyzed data span. The energy response and ef- 
fective area calibration are currently being validated, but are 
at present known to ~ 5% near 1 GeV, with uncertainties in- 
creasing to ~ 20% at < 0.1 GeV and ~ 30% at > 10 GeV. 

These data also substantially improve our knowledge of the 
pulse properties of the Vela pulsar (PSR B0833— 45) and are 
now placing important constraints on theoretical models. Al- 
ready in these data we see: 

1. The pulse profile is complex, with PI and P2 dominated 
by very narrow components and substantial structure in the 
'bridge' region. 

2. Although the PI and P2 phases are very stable across the 
7-ray band, the the P1/P2 ratio decreases with energy. There 
is a distinct third peak in the 'bridge' component which sharp- 
ens and moves to later pulsar phase with increasing energy. 

3. While faint emission appears within ~ 1° of the pulsar 
in the off -pulse phase, association if any with Vela is still un- 
clear, due to the bright Galactic background in the vicinity. 
We thus quote an upper limit on the unpulsed point source 
flux, at 2.8% of the phase-averaged pulsed flux. A true PWN 
association will require additional data, allowing resolution at 
higher energies and a match to the radio and TeV images. 

4. The phase-averaged 7-ray energy spectrum can be 
represented by a power law, with a exponential cut-off at 
2.9 ± 0.1 GeV. The hyper-exponential cut-off index b = 
0.88 ± OMtoit 

is not significantly different from the sim- 
ple exponential value b = 1. Large values of b, as expected 
for models radiating from the near-surface polar cap zone, are 
excluded. 

During the first-year sky survey, we expect to coUect ap- 
proximately 1.3 X 10^ additional pulsed Vela photons. This 
will allow a substantial increase in sensitivity to narrow pulse 
components. We also expect improvements in our under- 
standing of the instrument that will allow us to significantly 
reduce systematic errors. These two improvements wiU al- 
low derivation of high accuracy spectral properties in more 
than 100 bins of pulsar phase. These are the observational 
results that should let us pin down the location and energy 
distribution of the radiating particles. In turn this should help 
us 'reverse-engineer' this 7-ray machine and should provide 
substantial insight into the physics of pulsar magnetospheres. 
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TABLE 2 

Phase Averaged Spectral Points for the Vela Pulsar 



Energy (GeV) 


Differential flux F (erg cm •^s ^ ) 


AFstat 




AFtot" 


0.14 


1.381E-09 


5.805E-11 


2.952E-10 


3.009E-10 


0.18 


1.385E-09 


5.191E-11 


1.601E-10 


1.683E-10 


0.22 


1.514E-09 


3.895E-11 


2.261E-10 


2.294E-10 


0.28 


1.734E-09 


4.212E-11 


2.301E-10 


2.340E-10 


0.35 


1.953E-09 


4.563E-11 


2.232E-10 


2.279E-10 


0.45 


2.129E-09 


5.066E-11 


2.625E-10 


2.674E-10 


0.56 


2.088E-09 


5.289E-11 


1.671E-10 


1.752E-10 


0.71 


2.200E-09 


5.711E-11 


2.755E-10 


2.814E-10 


0.89 


2.405E-09 


6.505E-11 


3.006E-10 


3.076E-10 


1.12 


2.417E-09 


7.111E-11 


2.185E-10 


2.298E-10 


1.41 


2.493E-09 


8.007E-11 


4.626E-10 


4.695E-10 


1.78 


2.312E-09 


8.507E-11 


4.101E-10 


4.188E-10 


2.24 


2.273E-09 


9.770E-11 


4.164E-10 


4.277E-10 


2.82 


2.005E-09 


1.055E-10 


3.896E-10 


4.036E-10 


3.55 


1.617E-09 


1.052E-10 


4.524E-10 


4.645E-10 


4.47 


1.357E-09 


1.135E-10 


4.286E-10 


4.434E-10 


5.62 


1.019E-09 


1.116E-10 


1.854E-10 


2.164E-10 


7.08 


4.403E-10 


8.276E-11 


1.931E-10 


2.101E-10 


8.91 


4.435E-10 


9.914E-11 


1.870E-10 


2.116E-10 


11.22 


9.39()E-11 


6.150E-11 


O.OOOE+OO 


2.542E-11 


14.13 


4.191E-11 


4.550E-11 


O.OOOE+00 


3.043E-11 


17.78 


1.169E-10 


8.417E-11 


1.053E-11 


8.483E-11 


22.39 


O.OOOE+00 


O.OOOE+00 


O.OOOE+00 


O.OOOE+OO 


28.18 


O.OOOE+00 


O.OOOE+00 


O.OOOE+OO 


O.OOOE+OO 



"AFsyst is the systematic error on the flux and AFsj^^t are the statistical errors on the flux, while AFtot is the total error on the flux 



Data table for Electronic version. 



